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We study the structure of vacuum entanglement for two complimentary segments of a linear 
harmonic chain, applying the modewise decomposition of entangled gaussian states discussed in [23 |. 
We find that the resulting entangled mode shape hierarchy shows a distinctive layered structure 
with well defined relations between the depth of the modes, their characteristic wavelength, and 
their entanglement contribution. We re-derive in the strong coupling (diverging correlation length) 
regime, the logarithmic dependence of entanglement on the segment size predicted by conformal 
field theory for the boson universality class, and discuss its relation with the mode structure. We 
conjecture that the persistence of vacuum entanglement between arbitrarily separated finite size 
regions is connected with the localization of the highest frequency innermost modes. 



I. INTRODUCTION 

The study of entanglement properties of a number of 
physical models including spin chains, coupled fermions 
and harmonic oscillators [1-10], has revealed interesting 
aspects of entanglement in spatially extended many body 
systems. 

On one hand, a number of universal aspects connected 
to the behavior of the two-point correlation function has 
been verified. For instance, it has been shown that one- 
dimensional XY and Heisenberg spin chains near a criti- 
cal point Q , leads to the same scaling behavior for mass- 
less boson and fermion universality classes as predicted 
by conformal field theory [H El 13. In the massless 
case, the entanglement between a region of size L and 
the remainder of the system, increases either as ~ i log L 
(bosons) or ^ ^logL (fermions). These characteristic 
behaviors relate to a one-dimensional version of the black 
hole entropy "area-law" |T3. ITsllTfil l . 

In contrast to these universal results, other aspects of 
many body entanglement have proven to be model depen- 
dent and not entirely captured by the behavior of the cor- 
relation functions. For example, the entanglement length 
defined in ^3 as an analog of the correlation length, has 
been shown for a family of models be infinite while the 
correlation length is finite [isL IT^ . The converse situa- 
tion has also been demonstrated for spin chains P, 0| and 
harmonic chains '2^. In these cases, while the correlation 
length can be infinite or large, the entanglement between 
two sites truncates to zero for a relatively small separa- 
tion. It is therefore fair to say that we still lack a generic 
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characterization of entanglement in many body systems. 

In this paper we concentrate on an aspect that has re- 
ceived relatively little attention, but which could poten- 
tially sharpen the emerging picture. This is the entan- 
glement structure dictated by many-body interactions- 
the connection between the form of the Hamiltonian, the 
quantum state entanglement structure, and the spacial 
distributions associated with the quantum states. Our 
aim is to provide a detailed analysis of the bipartite 
ground state entanglement structure of the linear har- 
monic chain. 

In general the determination of the entanglement 
structure of a given quantum state is a complicated prob- 
lem; however, things are greatly simplified for the vac- 
uum state of an oscillator chain, which is a pure state 
of the Gaussian family |2lj |. An important feature of 
multi-mode pure Gaussian states is their fundamentally 
simple structure with respect to bipartite entanglement: 
it can be shown that pure Gaussian state entanglement is 
equivalent to products of entangled pairs of single modes 
[22. , so that the total entanglement is the sum of the 
1x1 modewise entanglement contributions; therefore, 
the canonical structure of gaussian bipartite entangle- 
ment is 1 X 1-mode gaussian entanglement. We therefore 
aim to characterize the spatial structure of the entangled 
modes within each of two complementary regions of the 
harmonic chain and connect this spacial structure to the 
corresponding entanglement contributions. This analy- 
sis is performed for both the weak and strong coupling 
regimes. 

An extrapolation of our results to the continuum limit 
shows agreement with previously known results, such as 
the |lnL entanglement behavior for bosons, as well as 
provides new insight into the entanglement characteris- 
tics of the vacuum. In particular, it shows that the in- 
clusion of an ultraviolet cutoff-in any way necessary to 
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regulate the corresponding massless relativistic quantum 
field theory in the presence of interactions-gives rise to 
a localization of the highest frequency modes around the 
midpoint of the block. Thus, it appears that contrary to 
the behavior of correlations, the long-distance behavior 
of entanglement between localized regions is directly con- 
nected with the high frequency modes. This somewhat 
paradoxical feature of field entanglement was noticed ear- 
lier in connection with the extraction of entanglement 
from the vacuu m [24| and violation of Bell's inequalities 
in the vacuum [23 . 

The article is organized as follows. In Section II. we 
begin with a short review of Gaussian states and the 
mode- wise structure of entanglement. We then study the 
structure of the local mode shapes within a block and 
their form by introducing a mode "participation func- 
tion". In Section III. we specialize to the case of the 
linear harmonic chain ground state, and study in some 
detail the dependence of the correlation functions on the 
coupling strength and size of the chain. We show that 
three regimes of behavior can be identified for a chain of 
fixed size and varying coupling strength. In Section IV. 
we present a simple application of the mode mapping for 
the case of a single oscillator with respect to the rest 
of the chain, and demonstrate how in this simple case, 
the properties of the correlation function show up in the 
behavior of the entanglement entropy. Section V surveys 
the main results of this paper on a qualitative level based 
on a numerical analysis of the modewise structure in the 
weak and strong coupling limits. This survey is comple- 
mented in section VI with an analytic study of the spatial 
modewise structure as well as the derivation of the ^ In L 
"area" law. Our results are compared with previous re- 
lated work in Section VII. Finally we conclude in Section 
VIII. In the Appendix we provide a connection between 
the discrete and continuum correlation functions. 



II. GAUSSIAN STATE MODEWISE 
ENTANGLEMENT 

The ground state of a linear chain of oscillators is a 
pure Gaussian state . Gaussian entanglement is char- 
acterized by the following simplifying property (see 
also [2^ for an alternative proof): If \iI>)ab is any Gaus- 
sian pure state of N modes entangling modes in two re- 
gions A and modes in B, then lip) may always be written 
as a product of two-mode and one-mode states 

Wab = 1^1),^,^ m,^,^ . . . 1^.),^,^ |0),^ |0),^ (2.1) 

where the IV'i)^ ^ are entangled states of one mode from 
the set A and one mode from the set B and 10) ~ and 
|0)g are products of vacuum states for the remaining 
modes. The fact that any pure Gaussian state can be de- 
composed according to H2.1|l implies that the bi-partite 
entanglement of the state is the sum of the entanglements 



from each one of the participating pairs. In turn, the en- 
tanglement of each pair is the von-Neumann entropy of 
the reduced density matrix obtained from the pairwise 
state \ipi) ^ ^ . In addition to the quantification of the 
total amount of entanglement, it also becomes relevant 
to investigate the relative contribution of the individual 
entangled modes, together with the relationship that may 
exist between the "shape" of these modes and their en- 
tanglement contribution. We use this section therefore 
to review and introduce some general technical aspects 
of pure entangled Gaussian state analysis. 



A. Review of Gaussian States 

We begin by reviewing some basic facts about Gaussian 
(pure or mixed) states. Let us represent the canonical 
variables of an A'^-mode system by the vector 

V^{q,pf, (2.2) 

where q = (gi, 92, gw)^, P = {pi,P2, ■■■,PnY ■ The 
commutation relations may thus be expressed as 

[Va,rii3\= iJap (2.3) 
where J is the so-called symplectic matrix 

A Gaussian quantum state p for a set of N modes is 
uniquely characterized by the first and second moments 
of 77. In dealing with entanglement aspects of Gaussian 
states, a shift in the expectation value of the canonical 
variables corresponds to a local operation. Thus, it may 
be assumed throughout that that {?]) — 0. In such a case, 
the state is entirely characterized by the matrix of second 
moments, the so-called phase-space 27V x 2N covariance 
matrix (CM): 

M = Rc(7??/^). (2.5) 

Of particular interest is the group of transformations 
preserving the Gaussian character of the state. Within 
the family of states with (77) = 0, the group is the ho- 
mogeneous group of linear symplectic transformations 
S £ Sp{2N, R) preserving the commutation relations un- 
der 77 = St], or equivalently, preserving the symplectic 
matrix under the similarity transformation SJS'^ — J . 
Under a symplectic transformation, a Gaussian state 
characterized by a covariance matrix M gets mapped to 
a Gaussian state characterized by the covariance matrix 
M = SMS'^. 

Somewhat analogous to the construction of normal 
modes for a linear system is the construction of modes 
in which the covariance matrix takes a particularly sim- 
ple form. A theorem due to Williamson |2^ states 
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that a certain symplectic transformation Sw always ex- 
ists that brings M to the normal form ( "Williamson nor- 
mal form" ) 

W = SwMS^^diag{Xi,X2,...XN,Xi,X2,-XN), (2.6) 

where the diagonal elements Xi, X2, ■■■Xn are referred 
to as the symplectic eigenvalues and must be greater 
than or equal to 1/2 according to the uncertainty princi- 
ple. The transformation Sw defines a new set of modes 
("Williamson modes") rj{m), with corresponding annihi- 
lation operators Um = '^'"y^^'" ■ terms of these modes, 
the state p may be written as a product of oscillator 
"thermal" states HI 



)(1 



(2.7) 



where hm = amflm is the Williamson number operator 
associated with the creation/ annihilation operators. The 
average number operator obeys a Bose-Einstein distribu- 
tion 
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2/3,„ _ 1 



(2.8) 



Since the covariance matrix is now diagonal with 
((9m)^) = {{VmY) = {nm) + ^, the symplcctic eigenvalues 
can be related to the average number operator according 
to 



X,n. — 



(2.9) 



and connected with the Boltzmann factor and /3,„ by the 
relations 



/3m = In 



A„, + 1/2 

Xm ~ 1/2 ' 



Am = icoth 



(2.10) 



(2.11) 



In the case of X^ = 1/2, the thermal state reduces to a 
projector onto the vacuum state |0)m annihilated by the 
destruction operator a^i ■ Otherwise, the state is a mixed 
state, with a von Neumann entropy S{Xm), where (using 
natural units) 

5(A) = (A+l/2)ln(A + l/2)- (A-l/2)ln(A-l/2) . 

(2.12) 



under local symplectic transformations bring the local co- 
variance matrices = Re{riAVA) ^^'^ = Rs{vbVb) 
into Williamson normal form. The decomposition is 
related with a general property of covariance matrices 
which satisfy the "isotropic condition" 



Xit. 



(2.13) 



a condition that is satisfied by the covariance matrix of 
any pure Gaussian state with Ao = 1/2. Partitioning the 
vector of all the modes SuS r] — r]A®VB, the full covariance 
matrix of the pure state may be written in block form as 



Ma K 



(2.14) 



where K — {r]Ai]B)- It is then possible to show |22| that 
as a consequence of the isotropic condition. Ma and Mb 
share the same symplectic spectrum in the respective 
sectors where the symplectic eigenvalues are larger than 
1/2. By performing local symplectic transformations 
VA — SaVa and r]B = SbVb bringing AIa and Mb to 
Williamson normal form, it is then possible to show that 
the obtained correlation matrix K = SaKSJ^ connects 
only the sectors in A and B with the same symplectic 
eigenvalue and vanishes on the elements with symplectic 
eigenvalue 1/2. This means that if the local symplectic 
spectrum is not degenerate (apart from the sector with 
symplectic eigenvalue 1/2), the transformed correlation 
matrix K connects only those pairs of modes in A and B 
with the same local symplectic eigenvalue. On the other 
hand, if there are degeneracies in the local symplectic 
spectrum, one can still perform an additional one-sided 
local orthogonal symplectic transformation that brings K 
to a form connecting the degenerate modes in a pairwise 
fashion. 

In this paper, we will only be dealing with Gaussian 
states for which the correlations between the q's and the 
p's vanish. This affords a particular simplification in the 
investigation of the entangled mode structure. For this 
we introduce the following notation for the coordinate 
and momentum covariance matrices: 



G = 



m 



H = 



m 



(2.15) 



In the absence oi q — p correlations, the local covariance 
matrices may thus be written in block diagonal form as 



Ma = Re{r]A Va) 



Ga 






Ha 



(2.16) 



B. Modewise Entanglement 

Turning now to the entanglement of pure Gaussian 
states, suppose the N modes are partitioned into two sets 
J] A and riB- Then, the particular set of modes in terms of 
which the decomposition I^TI is achieved are those which. 



The local symplectic spectrum can then be obtained from 
the square root of the doubly degenerate spectrum of the 
matrix — (JaMa)'^ (where Ja is the symplectic matrix of 
the A-modes), a matrix which is 



-(JaMa)' = ( ^-^0^^ Q^jj^ 



(2.17) 



4 



Thus, the symplectic eigenvalues are given by the square 
root of the eigenvalues of HaGa (or GaHa)- Therefore 
we can express a certain two-mode state in the decom- 
position (|2.1|l as 

n 

(2.18) 

The total bi-partite entanglement is then given by the 
sum of two- mode contributions ^i^m)- 

C. Mode Shapes 

The number states \n}^ in the above two-mode Gaus- 
sian state, are eigenstates of the Williamson modes num- 
ber operator Na^, which in turn can be expressed as 
a combination of the local modes tja- The question is 
then how do the local modes contribute to each of the 
collective Williamson modes? This can be answered by 
studying the symplectic transformation r/A = Satia be- 
tween the local and global modes. Particularly, we now 
wish to directly relate the symplectic transformation to 
the eigenvectors of HaGa and GaHa- First note that in 
the absence oi q — p correlations, the symplectic trans- 
formation will not mix the g's and p's, so we may write 
it as 

qa^Xqa PA=YpA XY^ = t (2.19) 

where the last condition guarantees that the transfor- 
mation is symplectic. Now define A to be the diagonal 
matrix with the symplectic eigenvalues of Ma- Thus A 
may be written as 

A = XGaX^ = YHaY^ . (2.20) 

Now, let u^™) and w'^™-' be right eigenvectors of HaGa 
and GaHa respectively corresponding to the symplectic 
eigenvalue (the m-th entry in A) A„i, so that 

{HaGaW""^ = A^7.(™) (G^i7^)i;(™) = A^„i;(™) , 

(2.21) 

and for simplicity introduce a normalization so that 

(u(") .t;(™)) = 1. (2.22) 

It is then possible to see that 

= CmGAU^""^ = d^HAv'-"'^ . (2.23) 

where the proportionality constants c^, must satisfy 
the condition 

Cmdm = . (2.24) 

A convenient choice to make is Cm — dm — A^^ , implying 
that M^™' and w^™) must be normalized so that 

(z;(™) . Hav^""'^) = («("' • Gau'"^) = Xra ■ (2.25) 



Next, we note that since Ga and Ha are symmetric, m^™-' 
and w'-™^ are respectively the right and left eigenvectors 
of HaGa respectively. This implies, together with our 
normalization convention, the orthogonality condition 

(u(™).^;(")) = ^„„ , (2.26) 
as well as the spectral decomposition 

if^GA = ^A^u(™)w(")^. (2.27) 

'm 

Using ifT^ and ^(TI^ . we finally arrive at the condi- 
tion that 

(i;(™) . HaV^'^'^) - («(") • G^?/(")) = \m5mn ■ (2.28) 

It is now straightforward to set up the symplectic trans- 
formation matrix X: Letting (jji be a column vector with 
all entries set to zero except the z-th one, one verifies by 
direct substitution that the matrices 

X = Y = Y,cj^raV^^^^ , (2.29) 

m m 

indeed satisfy 12.2()|l . Finally, we can now express the re- 
lation between the local and collective phase-space modes 
as 

D. Mode Participation Function 

Generically, the transformation connecting the initial 
modes to the Williamson normal modes is not an or- 
thogonal transformation. This means that in general the 
mode functions m*^™^ and w'^™) for the g's and the p's may 
be very different in shape. (For instance, it could be the 
case that the new g's may be fairly "localized" , with sig- 
nificant amplitude contributions from only a small num- 
ber of the old g's, while the new p's may show more a 
"collective" shape with more or less equal contributions 
from the old p's). Thus, it becomes convenient to de- 
fine a function taking into account, on an equal footing, 
the contribution from both the old g's and g's in a given 
Williamson mode. 

To this end, we note the expansion of a given 
Williamson mode creation operator in terms of local 
mode creation and annihilation operators. If the local 
site creation operator is defined as 

a^^\[^^q^+^7^P^] (2.31) 

where ^ is some arbitrary dimensional parameter, the 
creation operator for a given Williamson mode m on one 
side can be expanded as 

«™ - \ E(c^^f"' + e.^'fV. + (c^"!"^ - i^vt^M , 

(2.32) 
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In turn, this expression can be recast in the form 



i 

where 



cosh rj^™^ai + sinhr^^^al 



and 



-("^(i) = tanh" 



(m) , y-o (m) 



(2.33) 



(2.34) 



(2.35) 



We also note from the definition of the modes that 
. = 1 so that = 1- From the 

expansion of the mode creation operator, we see that the 
function V^™^{i) captures the weight of the of the local 
site participation, invariant under local rescalings and 
phase space rotations. We term it the mode participation 
function 



E. Mode Mapping 

We now wish to relate the the modes on the "A" -side 
with the modes on the "B" side corresponding to the 
same symplectic eigenvalue Xm > ^. This can be done 
by noting that from the isotropic condition (|2.13|) . the 
mode structure can not only be obtained from the local 
covariance matrices Ma or Mb, but also from the corre- 
lation matrix K. This in fact proves to be advantageous 
if the correlations between modes A and B are significant 
only for a small number of modes, as in the case of short 
range interactions. Assuming again no q/p correlations, 
the correlation matrix between A and B can be expressed 
in the block diagonal form 



K^ReivAvl) 



Gab 
Hab 



(2.36) 



where Gab = {qAqs) and Hab = (paPb)- Now, from 
(|2.13() one obtains the relations between Ma , M b and K 



{JaMaY ^{JaK){JbK^) 



1 



1 



{JbMbY + {JbK'){JaK) = -1 



(2.37) 



(2.38) 



where J a and Jb are the symplectic matrices in the A 
and B sectors respectively. Substituting in the forms 
()2.16(l and H2.36|l . we obtain the relations 



Km^ Gab 



^ H. 



AB 



Km^ Gab 



FIG. 1: Mode mapping. 



Next, defining nf^ to be the non-zero eigenvalues of 
—HabG^b^ thus see that the local symplectic eigen- 
values can also be expressed as 



A™ - ^ 



(2.41) 



Now, label the mode functions corresponding to the sym- 
plectic eigenvalue A„i > ^ as v^^^ and u^^\ w^-* 
for sides A and B respectively. According to H2.39|l and 
()2.21(l we see that they are solutions to the eigenvalue 
equations 



HabG, 



(m) 



lAB'^AB^A 

GabH^b^a^' 



and 



ilABGABUg 



(m) 



abJ^abVb 



(m) 



2 (m) 
1^771"^ A 

2 (m) 



.2 (™) 

2 (m) 
'^m'^B ■ 



(2.42) 

(2.43) 

(2.44) 
(2.45) 



Now, multiply both sides of say 12.42|l on the left by 



G^^, to find from lfn5|l that G^gHABiG^g 



u'T^), thus showing that w^"'' 



'T 

AB"-A 



Gab 



,M 



-K^„(G' 

A similar procedure applied to all the above equations 
shows that 



V^j^'> cx G 



■T , 
AB 



(m) 



CX H 



AB^A 



(m) 



(m) 
^A 
(m) 



CX Gabu 



zj (™) 

CX HABVg 



(2.46) 
(2.47) 



HaGa = 


1 

4 


- Hab Gab : 


(2.39) 


_ ^ f<T (rn) 
— ^AB^A 


(m 
^A 


HbGb — 


1 
4 


- HabGab ■ 


(2.40) 


,.(^^^) _ ^ }JT (m) 
Km 


(rn) 



The choice of the proportionality factors involved here is 
constrained by the normalization conditions imposed on 
the mode functions on both sides. A consistent choice 
preserving the normalization on both sides is 



^Gabu'^bH^AS) 
-A^Habv^bX'^AQ) 
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The above relations, together with the relations l|2.23|l , 
yield a systematic procedure by which all the mode func- 
tions on both sides, corresponding to a given symplectic 
eigenvalue can be constructed once a single mode func- 
tion is found. This proves particularly advantageous if 
the number of modes on one side is considerably smaller 
that in the other. This construction is summarized in 
figure 



III. HARMONIC CHAIN 

A linear harmonic chain of N local oscillators laid out 
on a circular topology may be modelled by canonical vari- 
ables {qi,Pi) with i = 1, .., TV, with the dynamics given by 
a Hamiltonian of the form 



H 



1=1 



p^ +<}i - a qniqi+i\ 



(3.1) 



where we identify Qn+i = qi and qo = gjv and the di- 
mensionless parameter a characterizing the strength of 
the coupling between adjacent neighbor sites is strictly 
less than unity. Note that such a Hamiltonian can be 
obtained from the standard Hamiltonian of a chain with 
spring-like nearest neighbor harmonic couplings 



1 ^ 

H = -y 

i=l 



M * 



(3.2) 



by means of the canonical variable re-scaling 



and the identification 



2K 



2K 



1 



2K 



2K 

(3.3) 



(3.4) 
(3.5) 



The last relation provides a restriction, < a < 1, to the 
possible values of the coupling constant in 1)3. The 
limit of strong coupling between neighboring oscillators, 
^j^2 ~^ oo, corresponds to a — > 1, and the weak coupling 
limit to a ^ 0. 

The Hamiltonian H3.1|l can be brought to a normal 
form by introducing a set of annihilation (creation) op- 
erators 0(6*^;) ( d^{9k) ) satisfying the commutation rela- 
tions 



(3.6) 



with the indexing angular variable 0k playing the role of 
a dimensionless wave number or pseudo-momentum and 
taking the values 

(^k = ^, (fc = 0,l,...,7V-l). (3.7) 



Defining the dispersion relation (in units of Eg) 



iy{9k) = Vl - a cos 9k 
and expressing g„ and p„ in the form 



(3.8) 



q-n 



Pn 



-y 



^M9k) 



v{9k) 



[a{,9k)t 



\a(Pk)e" 



■fh.c] , (3.9) 
h.c] , (3.10) 



the Hamiltonian (|3.1|) achieves the normal form 



H = EoY,^{ek) 



a\ek)a{ek) + 



(3.11) 



which is then diagonalized by the Fock states of the cre- 
ation/annihilation operators. 

In particular, we will be interested in the ground state 
|0), satisfying 



a(0fe)|O) -0, 



(3.12) 



for all 9k- For this state, the wave functions in the coor- 
dinate and momentum representations assume the Gaus- 
sian form 



■0o(g) oc exp 
V'o(p) oc exp 



^-^p^H-^p 



(3.13) 
(3.14) 



where the covariance matrices G and H , for Q and P 
respectively, satisfy the generalized uncertainty relation 
GH = i , with the entries defined by the two-point vac- 
uum correlation functions: 



Hi, 



(0|g,g,|0) = 5( 
(OIkPjIO) 



(3.15) 
(3.16) 



Furthermore, since the state is Gaussian, higher moments 
of the oscillator coordinates or momenta are expressible 
in terms of the two-point correlation functions. Thus, the 
relevant physical information associated with the vacuum 



is contained in the correlation functions 
which we now study. 



9(^3) 



and h 



A. Vacuum Correlation Functions and Three 
Regimes of Behavior 

The general entanglement behavior of the partitioned 
harmonic chain is dictated by the behavior of the cor- 
relation functions gi and hi defined in Eqs. H3.15() and 
1)3.16(1 . Their explicit form is given in terms of the dis- 
persion relation H3.8|l by 



■ cos{l9k) 



'l""^ ^ — ^K^fc)cos(?0/c), 



(3.17) 



(3.18) 
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and consistent with the translational invariance of the 
Hamiltonian, we note their dependence only on the sep- 
aration I = {i — j) mod N . In the hmit N ^ oo with 
a fixed, these expressions yield the Riemann sum for the 
integral of the argument as a function of a continuous 9 
ranging from to 2tt, with the replacement of the factor 
of in front by (27r)~^. The correspondence with the 
continuum one-dimensional scalar field theory can also 
be obtained from these expressions by taking appropri- 
ate limits, as discussed in the appendix. 

Now, for a fixed value of the strength parameter a and 
for sufficiently large values of the chain size N, the behav- 
ior of these functions in terms of I becomes independent 
of N and reproduces the N ^ oo behavior, which can be 
expressed exactly in terms of hypergeometric functions 
as 



(oo) 

9l - TT-, 



z 

2^1 



M2fJ^,1 + \,1 + 1,z'] (3.19) 



^(oo) ^ ^ 



flZ' 

2 



I 



2^1 



-i,Z-i,Z-fl,z2)(3.20) 



where 



are binomial coefhcients expressed in terms 
of Gamma functions and z and /i are given by 

1 



1- VT^^ 



(3.21) 



This behavior in the large N limit will serve as a reference 
in analyzing the behavior for finite N . 

The shape of the correlations gi and hi as a function 
of the separation I and the coupling strength is depicted 
in figures [3 and |3| for a fixed value of N . Similarly, we 
have plotted in^jthe behavior of their values at ^ = as a 
function of the coupling strength for different values of N . 
For these plots, we have found it convenient to introduce 
an ancillary hyperbolic angle ^, implicitly defined by its 
relation to the variables a and z: 



tanh^ a = tanh2^ 



(3.22) 



This becomes a convenient parameter as for small values 
of a, the appropriate expansion parameter is z itself, in 
which case z ~ ^ for ^ <C 1; similarly, as a approaches 
unity, ^ provides a logarithmic scale for this approach, 
with^~ iln(l-a)/2. 

Let us first analyze the case N — oo. For weak cou- 
pling, a ~ z <§; 1, the hypergeometric functions behave 
as 1 -f o(z^); thus, with the Stirling approximation for 
the binomial, we obtain the weak coupling behavior 



1 



m 



2^ 



1 



(3.23) 



for 1^1. The weak coupling correlation functions 
are therefore of short-ranged, exponential behavior, and 
characterized by the correlation length 
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FIG. 2: The vacuum gi — {qoQi) correlation function as a func- 
tion of / for different strength parameters where a — tanh 2^. 
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FIG. 3: The negative of the vacuum correlation function 
hi — {poPi) as a function of / (Z > 1) for different strength 
parameters. 



On the other hand, in the limit a ~ z — !■ 1 with fixed 
I, the gl°°^ correlation function diverges, and is deter- 
mined by the asymptotic behavior of the hypergeometric 
function 



Til + l) 



In 



1-z-^ 



+ V ( ) +7 



F(1/2)F(/ + 1/2) 

(l + 0(l-z)) 



(3.25) 



where 7 is the Euler-Mascheroni constant and ipin) is 
the DiGamma function. With the asymptotic expansion 
-0(/ + 1/2) ~ Inl for large > ^ > 1, this yields for 
z~l. 
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in the case of the correlation function hi, one finds the a 
independent limit 



lim hi = —'^——^ 

a^l TT 4^2 - 1 



(3.27) 
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FIG. 4: Variances in q and p as a function of coupling strength 
for different values of N. 



Next we consider finite-size effects. The first consid- 
eration has to do with the range the above limiting ex- 
pressions given a finite size of the chain. In this respect, 
the exponentially decaying behavior is not expected to 
hold once the correlation length becomes of the order 
of half the chain size; similarly, the strong coupling ex- 
pansion for gi is valid whenever the function is positive. 
For sufficiently large N, one should expect these condi- 
tions to be met at values of a which are close to unity, 
in which case 1 — z ~ ■y/2(l — a). Then the conditions 
Ic < N/2 and In ((1 - z)iV/4) < 0, lead essentially to the 
same transitional condition between short-range and long 
range correlation; the transition from short -ranged to 
long-ranged behavior then happens when a transitional 
correlation size scale, which we define as 



Nt{a) 



1 



(3.28) 



becomes of the order of N. Note incidentally that for 
Nt ^ I, the plotting parameter ^ is related to Nt ac- 
cording to 
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IniVt, 



(3.29) 



The second consideration has to do with the magnitude 
of the correlations for a close to unity, which is set by the 
the diverging qq correlation function and specifically its 
behavior at / = 0. For finite but large N, the value of gi 



can be approximated by the N = oo expression plus cor- 
rections in powers of arising from the error between 
the sum H3.17|l and the corresponding Riemann integral 
obtained when N oo. The relevant correction comes 
from the contribution of the 9k = mode in H3.17|l . 
which is not accounted for in the Riemann integral when 
the limit a —^ 1 is taken first. This contribution yields 
an 0{N~^) additional term, so that for a ^ 1, the finite 
size gi is approximated by 



{N) _ (oo) 
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(3.30) 



Since the correction arises from the contribution to the 
sum from the Ok = Q eigenmode, it accounts for a strong 
collective effect due the finite size of the chain. This 
correction dominates the coupling strength behavior of 
the gi correlation functions above a critical value of a, 
or, for fixed a, below a critical value of N . This criti- 
cal value is determined by requiring that the correction 
term should become of the magnitude of the logarithmic 
term depending on the coupling strength in the strong 
coupling N — oo correlation function, i.e. , 
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Taking a close to unity yields finally a critical chain size 
value Nr 



In TV, (a) 



(3.32) 



The above considerations allow us to distinguish three 
regimes of behavior for a chain size of a large and fixed 
value of N , for convenience to be referred to as the type 
/, // or III regimes: the type /, weak coupling regime, 
is determined by the condition A'' ^ Nt(a) > Nc{a), 
and is characterized by short-ranged correlations; when 
the condition Nt{a) ^ N is met, we enter the interme- 
diate, type // long-range regime, in which the scale of 
correlations is of order and In Nt (a) and the correlations 
decay logarithmically as a function of I; finally, the con- 
dition Nc{a) > N determines the long-range, type /// 
regime, with the same logarithmic decay as a function of 
the distance but in which the scale of the correlations be- 
haves as Nt{a)/N. The distinction between these three 
regimes will serve a guideline for the characterization of 
the different regimes of entanglement as well. 



IV. AN ILLUSTRATIVE EXACTLY SOLUBLE 
CASE 

Before proceeding with the numerical results in our pa- 
per, it may be useful to consider a simple example that 
is easy to solve and shows very general qualitative fea- 
tures of the dominant entanglement structure between 
two complementary regions of a chain. In this case we 
wish to understand the entanglement of one oscillator of 



the chain versus the rest of the chain. For such a parti- 
tion, it is quite easy to obtain the degree of entanglement 
by looking at the local covariance matrix of the single 
mode, which is given by 



Ma 



90 
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(4.1) 



The single symplectic eigenvalue of the matrix is A = 
\/ goho and thus the degree of entanglement between the 
oscillator and the remainder of the chain is given is given 
by S{^/goha) where S is defined by Eq. H2.12|l . 

To understand the behavior of the degree of entangle- 
ment we look at the relevant approximations for the three 
regimes discussed in section (|III All . First, expanding the 
hypergeometric functions in Ij3.19|l and (|3.2Q(I in a power 
series in z, we obtain for the symplectic eigenvalue the 
weak, type / regime. 
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FIG. 5: Entanglement of a single oscillator versus the rest of 
the chain as a function of coupling strength for a chain size 
— 10* (solid). The dashed lines show the approximations 
of Eq. 16. 2H : weak (I), strong transitional (II), and strong 
(III). 



Next, for strong coupling, Nt{a) > N, and assuming 
TV ^ 1, we approximate Hq and go by their limiting be- 
haviors for a —> 0. For ho the limiting value is V2/7r, 
while for go we can use the asymptotic ()3.25(l to show 
that 
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In (4 AT* (a)) , 



(4.3) 



where Nt{a) — ■\/2/(l — a) as defined previously. Thus 
the symplectic eigenvalue becomes 
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A// - -0n(47Vt(a)); 
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(4.4) 



this accounts for the type //, long-range regime with a 
logarithmic scale of the correlations. Finally, if, we take a 
large enough so that Nc{a) » N, then the go correlation 
function can be approximated by the single contribution 
of the collective 0^ = of the whole chain 
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1 fNt{a) 

2V2 V 



2NVT^ 2V2 V N 
This yields the symplectic eigenvalue 
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characterizing the strong long-range (type ///) regime. 

With these approximations we use the following expan- 
sions for the von-Neumann entropy 12.12|l as a function 
of A 



SiX) 



i«l 



(A - i) [1 - ln(A - i)] , A 

1 + In A , A > 1 



(4.7) 



This yields, for the characteristic behavior of the degree 
of entanglement in the three regimes, the expressions 



In^ 



1 + 1 In M2I 



Ntia) < N 



II : Nt{a) ^ N>Nc 
II : N^ia) > N 



(4.8) 

These three approximations compared quite with numer- 
ical computations of the entropy as a function of coupling 
strength in Fig for a chain size of TV = 10*. We note 
that the only importance dependence on the chain size for 
fixed coupling strength comes in the very strong regime, 
in which case the entanglement decreases logarithmically 
with TV. 

Correlated with the three regimes of behavior are the 
shapes of the u and v mode functions for the entangled 
mode on the complement of the chain. These are easily 
expressed in terms of the correlation functions gi , hi using 
the mode-mapping diagram Fig. ^ The l-th components 
are simply given by 



ui = 



\Jgoha 



\J90ho 



I = l,..,iV- 1, 



(4.9) 
(4.10) 



Thus, the shapes of u and v are easily inferred from fig- 
ures 121 and 121 in the weak regime, the mode functions 
exhibit rapid exponential decay away from the two sites 
adjacent to the singled-out oscillator; in the long-range 
transitional regime, the ui mode, being proportional to 
hi , decays as the inverse squared of the distance from the 
edge, while the u/ mode decays logarithmically; a similar 
behavior is exhibited in the strongest long-range regime, 
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FIG. 6: Total entanglement as a function of coupling strength 
(in terms of 5) and A^;, (log-2 scale) for a chain of size A*' — 
2048. 

except that the vi mode except that the logarithmic be- 
havior is essentially suppressed. Thus we see that in the 
strongest regime the entangled mode at the complement 
shows a collective correlated behavior in the oscillator 
momenta and a power-law correlation in the oscillator 
displacements. Two distinct limiting behaviors in a may 
then be inferred for the mode participation function of 
this mode, based on the asymptotics of the correlation 
functions, namely 

Pi-^l ■ (4.11) 



V. STRUCTURE OF GROUND STATE 
ENTANGLEMENT: QUALITATIVE STUDY AND 
NUMERICAL RESULTS 

We turn then to a qualitative study of the entangle- 
ment behavior of the harmonic chain vacuum and the 
underlying modewise structure, based on numerical re- 
sults. We shall be interested in a system of N oscillators, 
of which the first contiguous N}, < N/2 constitute a sub- 
system that we shall refer to as the "block", while the 
remaining N ~ Nf, will constitute what we shall refer to 
as the "complementary block". From the translational 
symmetry of the problem, the starting point of the block 
on the chain is completely arbitrary. 

Our first survey has to do with the overall behavior of 
the degree of entanglement of the block and its comple- 
ment as a function of the relevant properties: the size 
of the block Ni,, the size of the chain N, and the cou- 
pling strength. This overall behavior has been plotted 
in Fig. El for a chain of size N = 2048, sweeping the 
block size up to Ni, = N/2 and the coupling strength 
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FIG. 7: The three different qualitative regimes of entangle- 
ment as inferred from level curves of total entanglement in 
Fig. H 

so that Nc > N. As was seen in the previous section, 
the behavior of the degree of entanglement as a function 
of the coupling bears the signature of the three regimes 
outlined in section Ijlll A|) . This signature is also evident 
from Fig. [Hlwhen we look at slices of constant Ni,: As 
a function of ^ = tanh(2a), the entanglement rises from 
zero in the weak coupling regime to the characteristic 
plateau of the type // regime, which is more or less cen- 
tered around the value for which the transitional scale is 
of the order of the chain size (in our case around ^ ~ 3.2); 
as ^ is further increased, we see again the characteristic 
^ In Nt (i.e., linear in ^) behavior of the degree of entan- 
glement for the type /// regime. As we sweep then Nt 
and look at the entanglement as a function of ^, it then 
becomes possible to create an analog of a phase diagram 
on the £, — Nb plane for the qualitative behavior of the 
degree of entanglement; such a diagram is show in Fig. 
[3 and was constructed on the basis of the level curves of 
equal entanglement obtained to Fig. El 

Continuing with the behavior of the total entangle- 
ment as a function of the block size, we again find dis- 
tinct characteristic behaviors depending on the specific 
regime probed in the "phase diagram": for the type / 
regime, the total degree of entanglement remains essen- 
tially constant as Nj, is varied. Clearly, this is a manifes- 
tation of the short-ranged correlations characteristic of 
this regime, and suggests that entanglement between the 
block and its complement in this regime is due essentially 
to edge effects. More interesting, however, is the behav- 
ior in the strong type-// and type-/// regimes, where 
long range correlations are present and where therefore 
one should expect significant contributions from the bulk 
of the block. Indeed, one finds that in both of these 
regimes the degree of entanglement shows a logarithmic 
dependence on the size of the block-the one-dimensional 
analogue of the black hole area theorem. As shown in 
Fig. |S1 for sufficiently large values of the coupling and 
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FIG. 8: Entanglement as a function of InA^j, for A'' = 2048, 
5 = 12 fitted to a straight fine of slope 1/3. 
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InN, 

FIG. 9: Entanglement as a function of block size for different 
values of the chain size A*' ~ 3.5). 

N, and for <^ N the dependence takes the form 

E~Eo{a,N) + ^\nNb (5.1) 

reproducing the predictions of conformal field theory for 
the boson universahty class 0, 0| . It is also interesting 
to note the behavior of the value EQ{a, N) as a function 
of the chain size N for fixed a, which shows a similar be- 
havior to that of the entanglement of the single oscillator 
mentioned in the previous section. As shown in Fig. |51 
for large fixed value of a, the oc IniVf, curves decrease 
in height with N until a critical value of N; below this 
value the cx In Nh behavior saturates at a value that is 
independent of the chain size. 

We turn then to the question of "where does the en- 
tanglement comes from ?" , from the point of view of 
the modewise entanglement structure of the harmonic 
chain vacuum. In conformity with the different regimes 
of behavior entailed by the correlation functions gi and 
hi , different behaviors are reflected in the shapes of the 




FIG. 10: Site participation function for the entangled 
Williamson normal modes in a chain of = 32 -I- 48 con- 
tiguous oscillators, at four values of the coupling strength 
parameter (from top to bottom): a = 0.1, 0.6, 0.9, 1 — 10"''. 
The mode shapes are ordered front to back according to the 
symplectic eigenvalue of the respective mode, with the domi- 
nant mode at the front. 



Williamson modes, as we illustrate in Fig. ^| which 
shows the site participation function for all the entan- 
gled Williamson normal modes of both the block and its 
complement. Figure [TTl also shows a small sample of the 
mode shapes for the u^*") and w^™^ modes, for various 
values of m, and at two different coupling strengths. A 
number of general features may then be identified regard- 
ing the mode structure: 

The first feature, which is evident from Fig, ^] has 
to do with the parity of the mode functions u*^™^ and 
due to a parity selection rule to be discussed in 
section (jVI All . Indexing the modes as to = l,...iVb, 
in increasing order of their contribution to the total en- 
tanglement, we find that in all circumstances entangled 
modes of either the chain or its complement have a defi- 
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FIG. 11: A sample of Williamson mode shapes (solid for 
u''"'\ dashed for u'™') for a chain of size = 74 {N = 500), 
and for weak and strong coupling. Mode ordering is according 
to decreasing contribution to the total entanglement. 



nite parity (—1)™+^, under reflections with respect to the 
center of the respective block. Thus, for each entangled 
mode-pair, the modes on both sides of the chain have 
the same parity, and the mode-pair that contributes the 
most amount of entanglement always involves even mode 
functions on both sides of the chain. This feature shows 
that while it may be argued that the entanglement be- 
tween two regions of the chain is primarily due to edge 
effects, it nevertheless involves non-local behavior within 
each region. 

A second feature of the mode structure that is present 
in all circumstances is the existence of a characteristic 
turning point, henceforth denoted by which we de- 
fine as the location (measured from the center of the 
block) of the oscillator for which the participation func- 
tion 7']-™'' is maximal for the m-th mode. The general 
behavior of this turning point as a function of the mode 
number and coupling strength is shown in Fig. 1121 What 
is observed is that Z^'"-' is quite similar to the turning 
point of semiclassical solution to a wave equation: gen- 
erally, the mode participation shows an exponential-like 
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FIG. 12: Turning point location as a function of mode num- 
ber, for three values of the coupling strength (A'^i, = 64, 
N = 160). 



decay beyond the turning point, so that all the mode ac- 
tivity becomes effectively confined to site distances up to 
Zf") fr om the block midpoint; furthermore, insider this 
region, we find a less dramatic decay of the mode am- 
plitudes from the turning point, achieving the minimum 
amplitude at the midpoint. However, in contrast to semi- 
classical solutions, what we find is an anti-correlation 
between the amplitude and the local wavelength: the 
longest wavelengths are present at the midpoint, whereas 
the oscillations become of the order of the lattice spacing 
at the turning point. 

The last general feature that we shall discuss is ob- 
served more prominently as the coupling strength in- 
creases, and has to do with the typical wavelength of 
oscillation of the mode over the majority of the interior 
region bounded by the turning points. This characteris- 
tic wavelength decreases with the depth (m) from around 
the order of of the chain size, for the outer mode, to the 
lattice spacing for the innermost modes. Interestingly, 
we also find a "dual" structure to the mode shape oscil- 
lations if the u["^^ and Wj-'"'' mode entries are multiplied 
by the alternating factor (—1)* (corresponding to a shift 
of = TT in their discrete Fourier transform), and the 
modes are ordered from the inside out. Using the index 
convention n = Nh — m, we find that as the coupling 
strength increases, the mode shapes for the lowest values 
of n in this hierarchy of "demodulated" modes become 
virtually indistinguishable from harmonic oscillator wave 
functions of the corresponding value of n, with a variance 
of order VTVfc (Fig. I13|l . More generally, for the whole 
hierarchy the index n labels the number of nodes in the 
demodulated mode functions. 

With these considerations in mind, wc turn next to 



13 



n=Q y^^"^ 




n=1 J \ 

















FIG. 13: Mode shapes for the first three innermost modes 
for the strong regime case depicted in Fig. 1131 demodulated 
by the oscillating factor ( — 1)'. Here n = Nb — ra. 
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FIG. 14: Entanglement as a function of mode number in the 
weak regime, for different values of the coupling constant. 
Here TVj, = 12, TV = 500. 



the characterization of the regimes of behavior discussed 
earlier from the standpoint of the mode structure. Be- 
ginning with the weak regime, we find that the modes 
for both the block and the complement fall into pairs 
of even and odd combinations of oscillator sites that are 
precisely localized in increasing distance from the edge 
of the block (Fig. I12|l . so that a mode involving a cer- 
tain oscillator pair of the block is entangled with a mode 
of the complement involving the two oscillators' specu- 
lar images with respect to reflections about the interface. 
This "wedge-like" structure in the mode distribution of 
both sides of the chain observed in Fig. ^| can be un- 
derstood in the a limit, where the basis of localized 
oscillator sites serves an eigenbasis for the whole chain; 
then, the leading effect of any infinitesimal coupling be- 



tween contiguous oscillators is to recombine the localized 
modes so as to produce a mode basis consistent with the 
parity selection rules. We refer the reader to section lVI Al 
for more careful analysis discussion in this respect. 

As we increase the coupling strength, and hence the 
correlation length lc{a), the mode shapes become dis- 
torted to the same extent that the two participating sites 
lie within a distance from each other that is less than or 
equal to lc{a). Thus the distortion of the modes proceeds 
from the inside modes of the block and then outwards as 
the coupling strength is increased. This distortion in- 
volves a gradual diffusion of the mode participation to- 
wards the interior region bounded by the modes, as well 
as the establishment of a characteristic wavelength of os- 
cillation within this region. As the frequency spectrum 
of the modes becomes sharper, the width of the interior 
region of the mode becomes broadened, thus accounting 
for the gradual outward curving of the turning point lo- 
cation seen in Fig. ^1 

Together with this mode-shape distortion, we also find 
the behavior of the entanglement contribution of the 
modes to be dictated by the correlation length in the 
weak regime. In Fig. E| we have plotted on a loga- 
rithmic scale the dependence of the entanglement on the 
mode number in the weak regime. As stated earlier, in 
this regime the modes come in pairs of definite-parity 
combinations of oscillator sites at a precise distance from 
the edge; letting dm be this distance {dm = m/2 for m 
even, {m+l)/2 for m odd) we have found that the entan- 
glement is to leading order degenerate between modes of 
opposite parity and the same dm, independent of Nt (for 
Nf, > 2), and of exponential fall-off with dm according to 



Em(a) 



z{a) 



2(2d„,-l) 



2(2d„-l)ln(^) 



(5-2) 

where z{a) is as defined in eq. H3.21() : from the definition 
(|3.23|) of the correlation length Ic, the characteristic decay 
distance for the entanglement is then Zc/4. As the cou- 
pling strength is then increased, the degeneracy is lifted, 
with the innermost modes showing the greatest relative 
splitting. 

This "wave" of mode distortion and degeneracy lifting 
continues from the inside of the chain out as the coupling 
increases until the correlation length becomes compara- 
ble to the size of the block and the outermost modes 
are distorted. As the coupHng strength is increased be- 
yond this point, no appreciable change in the shape of the 
modes is detected. Thus a critical value of a determined 

by 



(5.3) 



sets a threshold beyond which the mode shape structure 
becomes frozen in its strong coupling configuration. It is 
this condition that underlies the transition between the 
type / and the type // transition in Fig. 

The onset of this transition and the ensuing behavior 
of entanglement thereafter is best appreciated from Fig. 
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1151 where we plot, as a function of the coupling, the total 
entanglement and the partial entanglement from the first 
four dominant modes. The main signature of the tran- 
sition is the lifting of the degeneracy involving the first 
(even) and second (odd) modes, which together up to 
that point are the predominant contributors to the total 
entanglement. Beyond that point, however, a clear de- 
coupling occurs between the first mode and the remaining 
modes: the first mode accounts for the behavior of the 
entanglement as a function of coupling, reproducing the 
three-regime curve (Fig. O observed in section (|1V|I for 
the case of the single entangled oscillator; on the other 
hand, the remaining modes become frozen in their cou- 
pling strength behavior. That this "freezing out" indeed 
occurs when lc{a) ~ Nb is best seen from Fig. which 
shows the level curves on the ^ — Nf, plane for the entan- 
glement of the first and second modes, together with a 
graph of the line lc{a) = iVf,. 

In this way we find that for strong coupling, the behav- 
ior of entanglement as a function of a, particularly with 
respect to the type // to type /// transition, is entirely 
due to the single outermost (to = 1) Williamson mode 
of the block. As observed from the mode shape profiles 
(Fig. Ill|l . this mode involves an essentially constant par- 
ticipation away from its turning points (located at the 
edge of the block) and may therefore be interpreted as 
a remnant of the 9k = collective normal eigenmode of 
the whole chain. In the context of the strong regimes, we 
shall therefore refer to this mode as the collective mode; 
similarly, we shall use the term residual modes for the 
remaining Williamson modes of the block. The reader is 

l|VI C|) respectively for 



referred to sections (jVIBII and 
further analysis of the collective and residual behaviors. 

Addressing finally the large block size and chain size 
behavior of entanglement, we find that in contrast to the 
large coupling behavior, the responsible modes for the 
~ IniVf, behavior are the residual modes. This can be 
seen from Fig. El which shows both the total entan- 
glement and the contribution from the first four domi- 
nant modes as a function of IniVf, for various coupling 
strengths. For large enough coupling, we find that the 
collective mode shows a "freezing" in its A'';, behavior, 
while the greatest variation with Nb is shown by the first 
couple residual modes. A more careful examination of 
the large Nb behavior, however, indicates that the en- 
tanglement contribution of the outermost modes grows 
slower than In A'b-rather like InlnA^f,. Thus, the InA^ 
behavior of the total entanglement has to be attributed 
to a cumulative effect from a certain number of residual 
modes that contribute significantly in the large A^;, limit, 
which turns out to be of order InA^b. If the modes are 
labelled by their characteristic frequencies, this phenom- 
ena translates to an enhancement by a factor In A";, of the 
density of states at zero frequency 29], and corresponds 
in the black-hole analogy to the divergence of the den- 
sity of states outside the horizon in the absence of a UV 
cutoff 

In this respect, it is instructive to look at the entan- 




FIG. 15: Contribution of the first four dominant modes (thin 
lines) to the total entanglement (thick line) as a function of 
coupling strength Nt for different values of the the block size 
(iV = 2048). 





FIG. 16: Level curves of entanglement contribution from first 
(left) and second (right) Williamson modes to the total entan- 
glement depicted in Figures |S| and [7| The dotted line shows 
the correlation length Ic, Eq. 13.2311 . as a function of coupling 
strength. 



glement contribution of all the residual modes, for which 
Nb < N/2 in the limit a ^ 1. As shown in Fig ((TH), the 
logarithm of the entanglement of the m-th mode shows 
in this limit a scaling behavior 



lni?„(Afc)~-Affc/( — 



(5.4) 



where f{x) is some nonlinear function proportional to 
the central frequency of oscillation of the mode. Now, 
the significant contribution to the entropy comes from 
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FIG. 17: Contribution of the first four dominant modes (thin 
hnes) to the total entanglement (thick line) as a function of 
A^i, for different values of the coupling strength (A'^ — 2048). 



modes for which In Em is of magnitude unity or smaller, 
that is, / < As we then show in section (IVIC|) . the 
function / behaves for small values of its argument like 



^/(a;) ln/(x) X . 



(5.5) 



This behavior implies that an outer layer of residual 
modes, with mode numbers m < InA^t, yield the rele- 
vant contributions to the entropy. 

To summarize the results of our qualitative survey of 
the bipartite entanglement of the harmonic chain, we 
emphasize again the three-regime framework depicted in 
Fig. [7| We observe a weak coupling regime (type /), 
characterized by short range correlations, well localized 
modes of definite parity, and degenerate entanglement 
contribution from even and odd combinations of oscilla- 
tor sites at the same distance from the edge. When the 
correlation length becomes of the order of the chain size, 
the mode shapes acquire a more collective behavior, a 
characteristic wavelength is established for each mode, 
and the degeneracy between contiguous modes of oppo- 
site parity is lifted. A distinction also emerges between 
the outermost, or "collective", mode and the remaining 
("residual") modes. The collective mode shows a min- 
imal contribution to the large chain size entanglement; 
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FIG. 18: Scaling curve of entanglement versus mode number 
for the strong regimes, for A*';, <^ N/2. The plot includes data 
taken from a chain of size A*' — 256 {Nt ranging from 4 to 32), 
and chains of size A'' = 1024 and A'' — 4096 (A^i, ranging from 
4 to 64), all using ^ = 10 {N^ ~ lO'^). 



however, it accounts for all the strong coupling constant 
behavior thereafter, including the transition between the 
type // and type /// regimes. Conversely, the residual 
modes become frozen in their large coupling entangle- 
ment contribution, but are responsible for the large chain 
behavior. This behavior is conjectured to come from a 
cumulative effect from a layer of the first ~ In Ni, residual 
modes. 



VI. MODE STRUCTURE: ANALYTIC RESULTS 

We complement in this section the qualitative survey 
of the previous section with a more careful treatment of 
some of the results presented there. Our main aim is to 
produce simple analytic models that capture the essential 
elements of the modewise entanglement structure in weak 
and strong coupling. 



A. Entanglement in Weak coupling 

As shown in our qualitative survey, entanglement for 
weak coupling involves modes that are definite parity su- 
perpositions of local oscillator-site modes. It will there- 
fore be convenient to briefly discuss the role of parity in 
our problem. 

Parity selection rules follow from the fact that while 
the local correlation matrices {Ga,Ha) and {Gb,Hb) 
lose the translational symmetry of the corresponding ma- 
trices for the whole chain, they still possess a symme- 
try with respect to reflections about the central indices, 
i.e., Gni,+i-i,Ni,+i-i = GNb+i-i,Nb+i-i- Thus, the G 
and H matrices can be written in block diagonal form as 
Ga = Ga"*" © Ga~ and Ha = Ha^ © Ha~, in terms of 
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blocks acting on subspaces of definite parity. Similarly, it 
is easy to show that with respect to the mode-mapping 
between the Williamson modes of the block and those 
of its its complement, it is also possible to show that the 
matrices Gab and Hab map modes of a given parity def- 
inite parity in the block to modes of the same parity in 
the complement. Thus, the whole problem of finding the 
Williamson modes in the size Ni, block, and their corre- 
lated counterparts in the complement can be reduced to 
two separate problems of finding Williamson modes for 
the sectors of size involving covariance matrices of 

definite parity. 

/,From Fig. E]it is evident that for small values of a 
the entanglement contribution decays exponentially with 
the mode number. This leading order exponential be- 
havior may be understood from simple arguments based 
on the "wedge" shape in Fig. ^| If for a given en- 
tangled mode-pair the (properly-normalized) mode func- 



tions M^'', 



and v'^jp^ are known, then the re- 



spective symplectic eigenvalue satisfies 



= vf~^^HAGAuT = 4 - vT' HabGbaUj, . 

(6.1) 

Using the resolution of the identity 1^ = X]/3 ''^^^'"b'''^ 
this yields 



(m) 



1 



(m)T , 



(m) 



E 



HabVb u'£' GbaUa 



(m) 



(6.2) 



However, we recall that Gba maps u"^^ to meaning 
that the only term surviving in the sum is the one in 
which f3 = m. Thus we have that 



{m)T 



„ (m) {m)Try 

HabVb Ub Gbau 



(m) 



(6.3) 



to correlations between a site in the block at a depth 
dm from the interface and a "mirror image" site in the 
complement at the same distance to the interface. 

A careful symbolic computation of this leading order 
behavior for small TVf, and assuming N — oo shows that 
in fact the leading-order behavior of the symplectic spec- 
trum is given by 

A™(z)-i = (J) +o{z''--'). (6.7) 

This leads to a mainly exponential dependence with log- 
arithmic corrections expressed in eq. 15. 2|) for the en- 
tanglement per mode in the weak coupling regime, and 
yield the solid lines in Fig. I14L showing excellent agree- 
ment for weak coupling up to values of z ~ 0.7. Our 
simple argument also leads us to expect from Eq. H6.6|l 
that the degeneracy between even and odd states in Eq. 
(16. 7|) should be lifted by an eigenvalue splitting of order 
^Nh+2rn~i^ thus bccoming more significant as the mode- 
depth increases, and with the even parity modes having 
the slightly larger symplectic eigenvalue. This qualitative 
behavior is in fact verified in Fig. 1141 



B. Strong coupling— The collective mode 

We turn next to the emergence of the collective mode 
in the strong coupling limit and its entanglement behav- 
ior. This has to do with the fact that the scale of the 
correlation function gi{a) diverges as a — *■ 1, together 
with the fact that the function behaves logarithmically 
as a function of I. Thus, the correlation function gi can 
be separated as 



9i = 5o(a) + A, 



(6.8) 



Now, in the weak coupling regime, we approximate the 
m-th even or odd mode functions in the weak regime 
by symmetric or antisymmetric combinations of localized 
site positions of depth dm, where dm = m/2 for m even 
and (m -|- l)/2 for m odd, as suggested by Fig. ^1 



where go is the diverging in a self-correlation function, 
and A/ is tends to fixed limit independent of a. In turn, 
the correaltion matrix Ga of the block takes the form 



GA^N^go xx^ + AG, 



(6.9) 



(m±) 

i 


(m±) 
= » 




_L X T fa A\ where AG 
i,d^ ± di^Nt-d„,+i\ (6.4) 


(m±) 


(m±) 
= ^B ^ 




'i,Nt+d„, ± Si^N~d,„ + l] (6.5) 



Then, using the symmetries of the correlation functions 
Eq. (ins yields 



A 



1 

4 



"-Ni ^ "'2d,„-l 



(N) , {N) 
dV ±92dL-l 



(6.6) 



Next, we make use of the weak limit expressions for the 
correlation functions, namely that both gi and hi behave 
as z', with hi negative for Z > 1. Since 2dm — I < Nj, 
always, then the leading order expansion of A — i is ex- 
pected to behave as z^i'^'^m-i) ^ Thus, in this approxi- 
mation, the symplectic eigenvalue is seen as being due 



1 



X 



1 1 



1 1 



Note that since for fixed Nf,, 

lim GaX Nbgoix + o{go^)] , 



(6.10) 



(6.11) 



the vector x becomes an eigenvector of Ga in that limit 
with eigenvalue Ntgo- 

Similarly, since the momentum correlation functions 
are regular and tend to a fixed value as a — s- 1 , the prod- 
uct HaGa is also split into diverging and finite parts: 



HaGa = m go{HAX){x^) + Ha^G . 



(6.12) 
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We now construct the vectors 

Uc cx Hax, "fc oc X, , (6.13) 
with normahzation set so that Uc ■ Vc ^ 1, and define 

= x^HaX = Nbgo , = x^HaX ■ (6.14) 



Then, we have that 

lim HaGaUc 
lim vJHaGa 



5x^xK + o(g/)] (6.15) 
gM^ + oig-')], (6.16) 



showing that Uc and Vc become the right and left eigen- 
vectors of GaHa in the strong coupUng Umit with sym- 
plectic eigenvalue 



Next, we compute the value of h^, 
1 



(6.17) 



EE 

j=i j=i 



TT 4(i-j)2 



671" 



^{Nb + -] + ln(4) + 7 



(6.18) 



where -0 is the DiGamma function. For large A^b, we 
approximate 

1 



bTT 



■ In (47Vfa) , 



so that 



A 



9o 
V2n 



In (4iVb) 



(6.19) 



(6.20) 



becomes the symplectic eigenvalue the approximation. 

The behavior of the collective mode entanglement for 
the type-// and type-/// regimes follows the same anal- 
ysis performed in section Hvl for the single entangled os- 
cillator, based on the two regimes of strong coupling be- 
havior of go determined by Nt{a) and Nc{a). Asymptot- 
ically, we find that 



E ' 



ilnlniVf, + ilnlnAft(a) , Nt{a) N > 
ilnln7Vfc + iln7Vt(a)/A^ , iVc(a) > iV 



(6.21) 

Thus, the entanglement curve of the collective mode is 
essentially that of the single entangled oscillator, except 
for sublogarithmnc corrections dependent on the chain 
size Nb. 



C. Strong Coupling— Residual Modes 

General qualitative and quantitative aspects of resid- 
ual mode entanglement in the strong regime can be il- 
lustrated through a simple analytical model in the con- 
tinuum, along the lines of similar models discussed in 



the context of geometric entropy in black hole physics 
hll and reduced density matrices for free electron chains 
[29l |. Such models are useful in deriving the correlation 
between mode number, wavelength and turning point lo- 
cation, and thus can account for the scaling relation de- 
picted in Fig 1181 and the density of states determining 
the InA^f, behavior. 

As discussed in section^ an eigenvalue problem equiv- 
alent to HGu — X^u for a given region A, is the eigen- 
value problem 



2 ('i 
K,„U A 



c 



HabGab 



(6.22) 



where Hab and Gab are the matrices containing corre- 
lations between the two complementary regions A and B, 
and where the eigenvalue Hm is related to the symplectic 
eigenvalue Am according to A^„ = j + k'^^ . We assume 
henceforth an infinite chain, and for the block adopt an 
index convention centered at the block midpoint, i.e., so 
that indices run from -(iVfc - 1 - l)/2 to (A^b - 1 - l)/2; 
no loss of generality is entailed by assuming iVf, odd so 
that indices are integers. With these assumptions, the 
matrix elements of T are given by the sums to infinity 
over both regions of the complement bordering the block 



hk-igk-] + E ^ 



k+igk+j 



(6.23) 



We base our computation of C on the asymptotic expres- 
sions for the correlation functions (|3.25|l and (|3.25l) , re- 
taining the terms yielding the leading order behavior in 
N(, upon converting the above sums to integrals. Thus 
we use 



1 



9^-J 



90 



V2Tr 
1 1 



2V2TT N- jf 



(6.24) 



where we include the diverging part go corresponding to 
the collective mode. Replacing the sums in (16.23(1 by in- 
tegrals in k from Ni,/2 to infinity, and defining the scaled 
variables x = i/Ni,, and y = i/Ni,, we finally obtain a 
result of the form 



Cii,j)=N-'rix,y)., 
where 



r r 



-F 



CA' 



(6.25) 



V2 



^ cs 



Ml - v^) 

X 



ln(iVi 



.90 



X2) 



47r2(i 



In 



1 + 2;/ 
l-2y 



.ln( ii±M(i:iMV (6.26) 



47r2(x-y) V(l-2a;)(l-l-2^)^ 
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The matrices Tcs and Tca are of rank one and deter- 
mine the behavior of the dominant mode for the even and 
odd sectors respectively; in turn, the bulk of the residual 
modes is determined by F^. 

To probe the solutions to the above eigenvalue prob- 
lem, we next consider the continuum generalization. 



1/2 



1/2 



dyT{x,y)ip{y) = -n^ilj{x) 



(6.27) 



For this, it will be convenient to perform a change of 
variables suggested by H6.31|l . mapping the interval 
(-1/2,1/2) to the real line 



= In 



1 + 2x 



V — \n 



1 + 2?/ 



1 - 2a; y V 1 - 2y 

and rescale the wave function according to 

■0(u) = icosh"^ ip{u) . 

The new eigenvalue equation then reads 



dv Tiu, v) ipi'^^) = — k^'0(m) , 



(6.28) 



(6.29) 



(6.30) 



where the new kernel T(u,v) is split in the form Tqs + 
FcA + Fi,;, with the respective terms given by 



V2 



i cs 



go 



V 

1 

Stt^ tanh 



■ tanh 



u - 



1 

u 
2 



In 



2 cosh (f) 



u , fu 
— tanh — 
2 V2 



V 2 ) 



(6.31) 



We now use F ~ F^ for the inner modes, which amounts 
to neglecting the lowest wavelength solutions. In this 
case, it is straightforward to see that the plane wave func- 
tion il}{u) = exp(ia;u) is an eigenfunction of F/?. The re- 
sulting integral for the eigenvalues can be performed by 
contour methods and yields 



'^'r — 



1 

8^ 



dv — r 

_oo tanh (I) 



:isinh"^(7rw) . (6.32) 



Finally, reverting to the variables cc, y and taking even/ 
odd combinations of the plane wave solutions, we find as 
finite- wavelength eigenfunctions of the continuum eigen- 
value problem H6.27|l . the solutions 



-4x2 1 sin(S'^(a;)) 



where 



Su{x) = a; In 



1 -|-2x 
1 - 2x 



(6.33) 



(6.34) 



The local wavenumber of these solutions, S'^ {x) , increases 
from the value 4a; as one moves away from the midpoint 
and diverges at the boundaries (x = ±1/2). Similar 
mode functions were obtained in '29| for the eigenmodes 
of the reduced density matrix of the free-electron chain. 
Note that the symplectic eigenvalues are 



1 



/I 1 -2 

A = W - + - sinh (ttuj) = - coth(7rcj) , (6.35) 

so the Boltzmann-like /3m factors associated with each 
Williamson mode according to eqs. H2.11|l are given by 
(3 — 2ttu!. For large f3, this yields an asymptotic expan- 
sion of the entropy E ^ e~'^'^'^, from which the scaling 
relation depicted in Fig. ITSl should then be expected to 
follow. 

We next connect these solutions to the corresponding 
solutions of the discrete chain. The first thing to note is 
that in the continuum approximation, the discrete modes 
correspond to averages of the continuum modes over the 
lattice spacing. Thus, the diverging oscillatory behavior 
of the continuum modes as the block edge is approached 
holds only up to a certain distance from the edge in the 
discrete approximation, corresponding to the point where 
the local oscillation wavenumber becomes of the order of 
the cutoff imposed by the lattice spacing. As one moves 
towards the edges beyond this point, the oscillations are 
washed out in the average and thus the mode amplitude 
decays. This is consistent with the fact, noted in sec- 
tion0 that at the turning points each Williamson mode 
shows oscillations of the maximum wavenumber. The 
turning points are therefore fixed by the condition that 



c 



dS_ 

dx 



(6.36) 



where C is number of order unity. With the definition 
(|6.34(l , this yields the relation between the turning point 
and u> 



Xt 



1-C 



4cj 



(6.37) 



Once the turning point is identified as a function of we 
can work out the quantization of the modes. For this we 
use the fact observed from our numerical calculations in 
section^ that when the interior modes are modulated by 
the oscillating factor (—1) '-corresponding to frequency 
shift by the cutoff frequency-the modes show a hierarchy 
in which each mode can be labelled by the number of 
nodes n + Nt — m. Writing the modulated functions like 



1 



1 - 4.x 



2 cos(7rA^f,a; - Suj{x)) , 
(6.38) 

and similarly for the antisymmetric mode, we demand 
that il^mod{x) show n nodes between the turning points; 
in this way we obtain the quantization condition 



TTNbXt{uj) - S^{xt{uj)) 



(6.39) 
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The resulting equation for to is best cast in a form that 
can easily be interpreted from the scaHng relation of the 
entropy H5.4|l . For this we define the dimensionless pa- 
rameters 



4w 



fi = 



(6.40) f 



where m = Nh — n labels the modes from the outside 
in, and / measures the central wavelength of the mode 
relative to the cutoff wavelength nNb; for large w, the 
entanglement of the mode is therefore given by 



(6.41) 



in accordance with equation H5.4|l . From the quanti- 
zation condition (|6.39|) . the function /(/z) is then the 
solution to the equation 



^ 2 Vi-VW7. 



(6.42) 



For the outermost modes, corresponding to small /i, small 
/, the relation becomes independent of Q to leading order 
and reduces to 



(6.43) 



as mentioned earlier. On the other hand, the fit for the 
innermost modes is sensitive to the precise value of ^. We 
have found that a value of C ~ 0.45 yields a remarkably 
good fit to the numerical data (Fig I19|l . as well as the 
turning point location depicted in Fig. ^1 

It is worth mentioning that in the present approach, 
the role of a cutoff frequency for the continuum is not 
only reflected in a cutoff of the resulting frequency spec- 
trum of the reduced density matrix, but also entails a lo- 
calization of the mode functions away from the interface 
between the block and the complement. The condition 
determining the turning point can only be approximately 
estimated from the continuum model and therefore fails 
to account for the exponential fall-off of the amplitude be- 
yond the turning point. The method also fails to account 



for the Aa 



-1/2 



scaled width of the the innermost 



mode functions , which is sensitive to the fall-off details. 



D. Residual mode contribution to the Entropy 

To obtain the asymptotic InA^f, behavior for Ni, <C iV, 
we note that the significant contributions to the entropy 
will come from modes for which lj is of order unity or 
smaller, since the entanglement is suppressed exponen- 
tially for cj > 1 ; since this implies that / will be of order 
1/Ni, or smaller, the approximation H6.43|l becomes more 
accurate for these modes as Nj, — s- oo. Thus, discard- 
ing terms of order unity in the logarithm, the eigenvalue 
equation for uj can be cast in the implicit form 



2 ln(A^b 



1 - 



(6.44) 



FIG. 19: Comparison between analytical and numerical re- 
sults for the scaled central wavelength / as a function of the 
scaled mode number, using C, — 0.45. 



Next, we can argue that as A'^^ — *■ oo, the term in brackets 
in the denominator of (|6.44l) can be neglected for the cal- 
culation of the entanglement; this term becomes impor- 



tant when TVi 



w <C 1 in the lower (since w > I, the 



approximation is already warranted in the upper end). 
For the lower end, we consider the smallest posible value 
of (jj, for m = 2; according to H6.44|l . this value is of order 
of 1 / In Nh and is therefore of an order N/, / In greater 
than the scale at which 16.44|) deviates from linearity. 
Hence we find the approximate linear expression for uj, 

(6.45) 



for the range of modes yielding relevant contributions to 
the entanglement as Nb —> oo. Using the approximation 
E{/3) ~ 1 — ln/3 for /3 <C 1, the entanglement of the modes 
for which IniVf, ^ m > 2, can then be estimated to be 



E„ 



In 



h^ 



(6.46) 



Consequently, the outermost residual modes yield a con- 
tribution of order InlnA'^t,, in the same way that the col- 
lective mode behaves for large Nb. 

To obtain the leading term in the asymptotic expan- 
sion of the the total residual mode entanglement, we note 
that for Nb ^ m \nNb the entanglement contribution 
is suppressed, the total residual entanglement can be es- 
timated to be 



where 



SW) 



dm 
'dp 



dfiS{l3) 



"'O 



5/3 _ 1 



- ln(l - e-0) , 



(6.47) 



(6.48) 



and where we evaluate the density of states assuming a 
linear relation aX [3 = lo = Q (note that this approxi- 
mation entails neglecting terms of order o(lnlnA'h) and 
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lower). With an integration by parts, the entropy can 
expressed as 

and from H6.45|l . the density of states is given by 
dm 1 dm 1 

— ^——^—\nNt. 6.50 

dp 2n duj TT^ 

The integral is standard in statistical mechanics and is 
given by 7r^/6. Thus we find that as Nj, — > oo, the resid- 
ual mode entropy is given by 

Er = UnNb +0(lnlniVb), (6.51) 

as expected. 

VII. COMPARISON WITH PREVIOUS 
RESULTS 

In this section we comment on the relation of our re- 
sults with earlier work. Audenaert et. al. J8| have stud- 
ied the entanglement in the circular linear chain model, 
for various choices of bipartite divisions. In this work, 
the logarithmic negativity Q has been used as a mea- 
sure for entanglement. Interestingly, we find that the 
von-Neumann entanglement seems to be, in the present 
problem, a more sensitive quantifier of the connection 
between entanglement and correlations. In Section Ilia 
we showed that the behavior of vacuum correlation func- 
tions quantifies three regimes: a weak coupling regime 
characterized by short-ranged correlations, an interme- 
diate regime that is reached when the correlation length 
is of the order of the whole chain, and finally, a long 
range correlation regime. The transition between these 
regimes is clearly manifested in the behavior of the von- 
Neumann entropy (see Figs. 5, 7, 15), but is absent in 
the behavior of the logarithmic negativity. The difference 
also shows up when comparing the dependence of the two 
measures on the total chain size N. In the particular case 
where Nf, = N/2 (Q, corollary 1), the logarithmic neg- 
ativity shows no N dependence. The von-Neumann en- 
tropy, on the other hand, decreases as a function of N like 
ln{Nt/N) (seeeqs. (4.8), (6.21) and Fig. 9). This depen- 
dence on iV, can be physically understood as due to the 
contribution of the collective mode to the entanglement, 
which reduces with increasing N. Thus is seems that the 
logarithmic negativity is not sensitive to the contribution 
of the collective mode. 

The entanglement of a finite region for a one dimen- 
sional field bosonic and fermionic fields has been previ- 
ously investigated in connection with the black hole en- 
tropy "area law" [T^ ITsl IT^ . By employing methods of 
conformal field theory it has been shown JT., J^J that 
in the massless case entanglement behaves like 



where e plays the role of the UV cutoff, and c and c are 
the holomorphic and anti-holomorphic central charges of 
the conformal field theory, with c = c = 1 for bosons, and 
c = c = 1/2 for fermions. Thus the overall coefficient is 
given by either 1/3 for bosons, or 1/6 for fermions. The 
same type of universal behavior has been recently derived 
the XY and Heisenberg spin-chain models 3]. 

In the present work, we obtained a logarithmic de- 
pendance of the entanglement ^IniV, corresponding as 
expected to a bosonic field (Fig. (jSl). Furthermore 
we have seen that the logarithmic increase of the can 
be understood as an increase in the number of relevant 
contributing modes while the coefficient 1/3 can be ob- 
tained from the density of the modes. As in previous 
results, these modes can be identified to be in a layer 
that becomes infinitesimally narrow in the limit of large 
N. However, no attention has so far been given to the 
structure of the inner modes. A central outcome of the 
present work is that the inclusion of an ultraviolet cut- 
off, which is needed for the consistency of the one di- 
mensional field theory, gives rise to a localization of the 
highest frequency modes around the midpoint of the re- 
gion. Although the contribution of these modes to the 
entanglement is exponentially small it is plausible that 
these inner modes play an important role in the persis- 
tence of vacuum entanglement between separated regions 
as we suggest in the next section. 

A number of results presented in this paper can also be 
related to previously obtained results for both fermions 
and bosons in the context of the densit y rn atrix renor- 
malization group (DMRG)!!!!!!!!!!!!]!!!!!. In par- 
ticular, the factorized thermal form of the reduced den- 
sity matrices, the shape of the corresponding mode func- 
tions in the continuum limit, and the approximately lin- 
ear behavior of the frequency spectrum for the outermost 
modes have been studied extensively in that context. It is 
possible that the method used in Section FVI CI whereby 
the cutoff imposed by the lattice spacing is used to estab- 
lish the turning points and quantization condition of the 
mode functions, may prove useful for the DMRG scheme. 



VIII. SUMMARY AND DISCUSSION 

Previous work on ground state entanglement in chain- 
like systems has mostly focused on the dependence of 
the amount of entanglement on parameters such as the 
block size, the separation between sites and the nature of 
the bipartite splitting. While in the present work we have 
reproduced several earlier results, as discussed in the pre- 
vious section, the emphasis here has been on the study 
of the spatial entanglement structure emerging from the 
modewise decomposition of the ground state wave func- 
tion. 

From this analysis, we have identified certain general 
properties of the mode structure and its relation to the 
entanglement contribution. A first feature is localization- 
a definite characteristic distance from the division inter- 
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face for the entangled modes at either side of the inter- 
face, thus estabUshing a characteristic distance separat- 
ing the entangled mode functions. This in turn, serves 
to characterize the strength of the entanglement, which 
decays exponentially with this distance. 

A second feature, which becomes sharper with increas- 
ing coupling strength, is a characteristic wavelength cor- 
related with the degree of localization of the modes. This 
correlation is in fact observed in two guises: on the one 
hand, it provides an alternative characterization of the 
modes in terms of their central wavelength, with the in- 
nermost modes possessing the shortest wavelength dic- 
tated by the lattice spacing; on the other hand, the am- 
plitude of each mode is correlated to the local oscillation 
wavelength, with the largest amplitude occurring when a 
the cutoff wavelength. 

We have shown that the effect of the interaction 
strength on the shape of the modes and their contribu- 
tion to entanglement is fundamentally connected with the 
correlation length. When the coupling is strong enough 
such as the correlation length becomes comparable to 
the size of the system, and the system becomes effec- 
tively massless, scale free behavior emerges for the bulk 
of the modes. The shape of the mode functions can be 
connected to the scale free continuum field theory, and 
both the localization and the characteristic wavelength 
scale with the size of the block. 

On a more speculative note, it is possible the re- 
sults of this paper may shed new light into several fea- 
tures of mixed state entanglement for separated non- 
complementary regions in vacuum. It has been shown 
that for arbitrarily separated regions, vacuum entan- 
glement persists and Bell inequalities are violated [2^ . 
with a lower bound of the entanglement that goes like 
exp(— L^/D^) where D and L denote the size of the re- 
gions and their separation. It was shown that a large 
probe energy gap is required in order to extract this en- 
tanglement. This seems to suggest that the localization 
of the inner modes and their short wavelength charac- 
teristic are linked with the persistence of vacuum entan- 
glement at large distances. It is possible this persistence 
represents an effective shielding of the entanglement con- 
tent of the innermost modes. This qualitative argument 
could help explain the truncation effect that takes place 
beyond a critical distance in the discrete version of this 
problem when the region sizes are kept fixed, and could 
explain the discrepancy between entanglement and cor- 
relation lengths in other models. 



IX. APPENDIX: DISCRETIZATION OF THE 
MASSIVE CONTINUUM THEORY 

In this appendix we connect our results for the dis- 
crete chain with a continuum bosonic theory in one di- 
mension. To this end, consider the Hamiltonian for the 
one-dimensional massive continuum theory on a circular 
topology 

H=^Jdx [nix)^ + {(p'{x)f + ^?4>{xf] , (9.1) 

where we assume that x runs from —L/2 to L and that 
(f>{x — L/2) = (j){x + L/2), and field configuration and mo- 
mentum operators satisfying the commutation relations 



{x,t),Tr{y,t)] = iS{x - y) . 



(9.2) 



The Hamiltonian is diagonalized in terms of normal 
mode creation (annihilation) operators a(fc„)(a(A;„)^) 
with fc„ — ^n, n G Z, and such that [a(A;„), a(fcm)^^] = 
<5n,m- The fields are then given by 

^{x) = -ifE Kfc«)e^''""+h-c.] ,(9.3) 



7r(x) 



[a(fc„)e''="^-h.c.] , (9.4) 



VI^ V 2 

where Lo{k) is the usual dispersion relation Lu{k) — 
y/i^ + A?. The continuum field correlation functions can 
then be obtained yielding 



~(L)f \ ^ cos(fe„a;) 



(9.5) 



2L ^-^ i-^ikn) 
U'^^x) = J-^c^(fc„)cos(fc„a;) (9.6) 



In the limit when L —i- oo, these expressions can be ex- 
pressed in terms of the modified Bessel functions: 



s(oo) 



ix) 



1 f°° „ C0s(fc2:) 
— / dk- 



47r 7_oo "'^ ^W+Jj^ 
Ztt 



(9.7) 



/i(°°)(x) = — / dky/WTJ^ cos{kx) (9.9) 



-^i^r(.N) 



(9.10) 
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The asymptotic form of the correlation functions is 
given for cc ^ /i^^ by 



1 



2tt 



1 



2tt\x\ 



(9.11) 
(9.12) 
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and for x ^ ^ ^ , we have 



-fj,\x\ 



2^27r^|j 



87r|a; 



(9.13) 
(9.14) 



The theory at the continuum may be approximated by 
a Unear chain of N sites in a ring topology, introducing 
discretized variables g„ andp„, which up to scale changes 
samples the field and its conjugate momentum field at 
points Xn = —L/2+jjL. To obtain a Hamiltonian of the 
form H3.1|) , we first approximate the field Hamiltonian by 
replacing / rfa; ^ ^ "^'(^^n) ^ '^'^^ 
we perform the following transformation 

0(x„) = ^ ^(2;„) = Ai/2p„ (9.15) 

It is can now be seen that with the choice 

1/2 



A = 



we obtain a Hamiltonian of the form H3.1|l with 
K = A 

and an effective coupling strength 



a{tiL/N) = 



1 



= 1 - 



(9.16) 



(9.17) 



(9.18) 



The correlation functions in the continuum are charac- 
terized by the length scale As can be seen from eqs. 
(19.1119.14(1 ,for separation x < l//i the theory behaves es- 
sentially as massless, while for x > 1//J,, the correlations 
decay exponentially. This length scale can then be re- 
lated to the correlation length obtained from the infinite 

harmonic chain by the relation jia 



We can then verify that the following relations hold be- 
tween the discrete and continuum correlation functions 



q^ix) = — = lim Q^^^ 



a{pL/N) 



h^ix) = V2 lim ( — 

JV^oo \ L 



"'n(x,N) 



(9.19) 
(9.20) 



In finite massive chains, another useful limits are 



\x) 



1 (oo) 

—pz lim q 



A/\/2 



/i°°(a;) = V2 lim A^ h^°°\ , ^ 



—= lim q] 

V2 i-oo ^' 



(oo) 



or equivalently 

}f^{x) = \/2 lim 



(oo) 



(9.21) 
(9.22) 

(9.23) 
(9.24) 
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